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Abstract 

Empirical  databases  of  crystal  structures  and  thermodynamic  properties  are  fundamental  tools  for  materials  research. 
Recent  rapid  proliferation  of  computational  data  on  materials  properties  presents  the  possibility  to  complement  and 
extend  the  databases  where  the  experimental  data  is  lacking  or  difficult  to  obtain.  Enhanced  repositories  that  inte¬ 
grate  both  computational  and  empirical  approaches  open  novel  opportunities  for  structure  discovery  and  optimization, 
including  uncovering  of  unsuspected  compounds,  metastable  structures  and  correlations  between  various  character¬ 
istics.  The  practical  realization  of  these  opportunities  depends  on  a  systematic  compilation  and  classification  of  the 
generated  data  in  addition  to  an  accessible  interface  for  the  materials  science  community.  In  this  paper  we  present 
an  extensive  repository,  aflowlib.org,  comprising  phase-diagrams,  electronic  structure  and  magnetic  properties,  gen¬ 
erated  by  the  high-throughput  framework  Aflow.  This  continuously  updated  compilation  currently  contains  over 
150,000  thermodynamic  entries  for  alloys,  covering  the  entire  composition  range  of  more  than  650  binary  systems, 
13,000  electronic  structures  analyses  of  inorganic  compounds,  and  50,000  entries  for  novel  potential  magnetic  and 
spintronics  systems.  The  repository  is  available  for  the  scientific  community  on  the  website  of  the  materials  research 
consortium,  aflowlib.org. 
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1.  Introduction 

The  development  of  new  materials  is  a  difficult  en¬ 
terprise  because  of  the  vast  number  of  possible  combi¬ 
nations  of  elements  and  crystal  structures.  In  develop¬ 
ing  a  new  material,  or  improving  one  already  in  use, 
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the  materials  scientist  usually  relies  on  the  huge  body  of 
empirical  information  on  crystal  structures  and  materi¬ 
als  properties  compiled  into  various  databases.  Among 
the  most  widely  used  databases  of  inorganic  materials 
are  those  of  metals,  intermetallics  and  minerals  crystal 
structures:  CRYSTMET  [1],  ICSD  [2],  Pearson’s  Crystal 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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Figure  1:  Input  interface  of  the  aflowlib.org  alloy  library  of  binary  systems. 


Data  [3],  Pauling  File  Inorganic  Materials  Database  [4]) 
and  binary  alloy  phase  diagrams  [5,  6].  These  empir¬ 
ical  databases  are  understandably  incomplete,  not  only 
due  to  the  vast  combination  space  (most  of  which  has 
not  been  explored  [7]),  but  also  since  experimentation 
in  many  cases  is  difficult,  requiring  high  temperatures 
and/or  pressures,  very  long  equilibration  processes,  or 
materials  that  are  hazardous,  highly  reactive,  poisonous 
or  radioactive.  Computational  compilation  of  materials 
properties  is  more  feasible,  and  will  lead  to  much  more 
complete  repositories. 

The  dramatic  advances  made  over  the  last  half  cen¬ 
tury  in  calculating  material  properties  quickly  and  ac¬ 
curately,  based  on  quantum-mechanical  approaches, 
have  spawned  the  new  discipline  of  computational 
materials  science.  They  make  it  possible  to  create 
computationally-derived  repositories,  that  significantly 
complement  existing  empirical  data,  and  use  them  for 


rational  material  design.  Generating  such  repositories 
requires  efficient  screening  of  large  sets  of  structures 
with  many  different  combinations  of  elements  and  com¬ 
positions  in  a  high-throughput  (HT)  fashion. 

In  this  paper  we  describe  aflowlib.org,  a  library  of 
databases  of  material  properties  generated  using  the  HT 
framework  Aflow  [8],  that  may  be  mined  a  posteriori 
for  various,  a  priori  unanticipated,  applications.  It  cur¬ 
rently  includes  repositories  of 

i.  phase-diagrams,  free  energies,  stable  and  metastable 
structures  of  binary  alloy  systems  (Binary  Alloy  Project , 
Section  2) , 

ii.  electronic  structure  of  inorganic  compounds  (Elec¬ 
tronic  Structure  Database ,  Section  3)  and 

iii.  electronic  and  magnetic  properties  of  alloys  (Mag¬ 
netic  Database ,  Section  5). 

The  databases  are  undergoing  continuous  development 
and  extension.  They  are  accessible  on  our  website, 


2 


# - 

NbRh/ 2  92 

#  Structure  PRE  - 

Nb_svRh_pv .292 
-66  . 076200 

3.99980000000000 

0 .00000000000000 

0 . 00000000000000 

0 .00000000000000 

3 . 99980000000000 

0 . 00000000000000 

0 .00000000000000 

0 .00000000000000 

5.06540000000000 

2  2  Direct  (4)  [A2B2] 

0.25000000000000 

0.25000000000000 

0 . 23420000000000 

Nb 

0.75000000000000 

0 . 75000000000000 

0 .76580000000000 

Nb 

0.25000000000000 

0 . 75000000000000 

0 . 00000000000000 

Rh 

0.75000000000000 

0.25000000000000 

0 . 00000000000000 

Rh 

#  Structure  POST  - 

Nb_s vRh_pv .292 

0 . 934225 

4.32163818537409 

0.00000000000000 

0.00000000000000 

0.00000000000000 

4.32163818537409 

-0.00000000000000 

0.00000000000000 

0.00000000000000 

4.15998906360566 

2  2  Direct  (4)  [A2B2] 

0.25000000000000 

0.25000000000000 

0 .49523191611840 

Nb 

0.75000000000000 

0.75000000000000 

0.50476808388160 

Nb 

0.25000000000000 

0.75000000000000 

0.00000000000000 

Rh 

0.75000000000000 

# - 

-0 . 43378830000000 

0.25000000000000 

0 . 00000000000000 

Rh 

#  Ef/at  (VASP) 

0 .00000000000000 

#  Mom/at 

15 . 83740000000000 

#  Volume/at 

0.50000000000000 

#  Ca 

0.50000000000000 

#  Cb 

P  4 / nmm  #12  9 

#  space  group  PRE 

P  4 / mmm  #123 
# - 

#  space  group  POST 

XX 

Figure  2:  Structure  data  example,  pre-  and  post-relaxation,  in¬ 
cluding  the  Aflow  structure  label,  unit  cell  vectors,  composi¬ 
tion  and  atom  coordinates.  In  this  example  NbRh-BlO  relaxes 
to  Ll0. 


aflowlib.org. 

Other  online  materials  initiatives  include  the 
Materials  Project,  (materialsproject.org)  [9],  and 
the  Computational  Materials  Repository  (CMR), 
(wiki.fysik.dtu.dk/cmr). 
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Figure  3:  Example  of  a  structure  list,  showing  NbRh 
equiatomic  structures  calculated  by  Aflow. 


[NbRh]  opar=675. 275757  calcs=228  date:  2011-09-12  [http://materials.duke.edu] 


2.  The  Binary  Alloy  Project 

An  extensive  database  of  binary  phase-diagrams, 
compounds,  and  metastable  structures  (the  “binary  alloy 
library”  project)  was  created  by  using  the  Aflow  frame¬ 
work  [8].  The  generation  of  the  library  required  several 
tens  of  millions  of  cpu  hours.  It  currently  contains  over 
150,000  alloy  entries  spanning  the  entire  composition 
range  of  more  than  650  binary  systems  of  the  transition 
metals,  metals,  alkaline  earths  and  alkali  metals.  Each 
entry  includes  the  groundstate  energy  and  magnetic  mo¬ 
ment  per  atom,  the  crystal  structure  space  group  and 
prototype  or  the  Strukturbericht  designation.  This  alloy 
database  is  updated  regularly  and  extended  periodically. 
It  may  be  accessed  online  at  aflowlib.org  through  the 
user  interface  shown  in  Figure  1 . 

The  interface  allows  the  user  to  review  the  following 
options  for  each  of  the  binary  systems  included  in  the 
repository: 

•  all;  The  complete  list  of  structures  calculated  for 
the  specified  system.  The  original  and  relaxed 


Figure  4:  Example  of  a  binary  phase  diagram,  NbRh,  gen¬ 
erated  automatically  by  Aflow  from  groundstate  energies  of 
over  220  crystal  structures. 

structure  data  (including  unit  cell  vectors,  atom 
coordinates,  calculated  energy,  magnetic  moment 
and  volume  per  atom)  is  provided  for  each  en¬ 
try,  as  shown  in  Fig.  2.  Fig.  3  is  an  example  of 
a  list  of  equiatomic  structures  calculated  for  the 
NbRh  binary  system,  where  Fl0  is  identified  as  a 
groundstate  of  the  system  (labeled  by  Emin-GND). 
It  demonstrates  the  ability  of  Aflow  to  identify 
relaxation  of  structures  into  others,  e.g.,  the  sec¬ 
ond  structure  listed  in  Fig.  3  starts  as  prototype 
B10  (space  group  #129)  but  relaxes  into  Flo  (space 
group  #123). 

•  bcc,  fee  or  hep;  lists  of  bcc-,  fee-  or  hep-based 
structures  calculated  for  the  selected  system,  re¬ 
spectively.  These  lattice  specific  lists  are  useful  for 
implementing  cluster  expansion  calculations  [10]. 
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Figure  5:  The  input  interface  for  the  crystallographic,  electronic  structure,  thermoelectric  and  magnetic  database  of  aflowlib.org. 


•  phase  diagram;  the  zero-temperature  phase  di¬ 
agram  of  the  selected  binary  system.  For  each 
structure  in  the  database,  the  formation  enthalpy 
was  calculated  from  the  lowest-energy  elemental 
form  of  the  constituents.  Aflow  automatically 
constructs  the  convex  hull,  from  the  structures  of 
extremal  energy  at  the  various  concentrations,  and 
generates  the  phase  diagram  as  shown  in  Fig.  4. 

This  database  of  binary  alloys  has  been  used  in  sev¬ 
eral  studies  to  identify  hitherto  unknown  stable  phases 
in  several  alloy  classes.  In  one  case,  34  Mg  bi¬ 
nary  systems  were  examined  [11].  In  about  one  third 
of  the  systems,  the  ab  initio  results  identified  all  the 
known  experimental  phases  and  predicted  additional 
stable  phases.  Moreover,  seven  systems  Mg-Na,  Mg- 
Tc,  Mg-Zr,  Ir-Mg,  Mg-Pb,  Mg-Pd,  and  Mg-Rh,  thought 
to  be  non-compound-forming  based  on  experimental 
evidence  were  predicted  to  have  stable  phases.  This  new 
thermodynamic  information  may  be  useful  for  the  de¬ 
velopment  of  improved  lightweight  alloys. 

High-throughput  studies  were  also  undertaken  for 
several  platinum  group  metals,  Rh,  Re,  and  Ru  [12, 
13,  14].  These  metals  are  of  particular  interest  because 
of  their  many  applications  as  catalytic  materials.  The 


picture  of  these  alloys  that  emerges  from  the  extensive 
calculated  database  is  considerably  different  from  that 
depicted  by  experimental  data.  Many  of  the  ostensi¬ 
bly  phase- separating  systems  are  predicted  to  have  sta¬ 
ble  structures.  Although  some  of  the  predicted  phases 
may  only  be  stable  at  temperatures  where  sluggish  ki¬ 
netics  inhibits  their  formation,  they  serve  as  a  guide 
to  further  experimental  studies,  usefully  narrowing  the 
field  of  compositions  to  be  attempted.  Furthermore, 
even  if  the  bulk  phases  might  not  exist  at  room  tem¬ 
perature,  ordered  structures  might  be  stable  at  the  nano¬ 
scale,  e.g.  for  a  few  nano-catalysts  (Fe-C,  FeMo-C,  Pt), 
it  was  shown  that  the  surface-tension  contribution  to  the 
free  energy  plays  a  fundamental  role  in  stabilizing  or¬ 
dered  structures,  thereby  causing  catalytic  deactivation 
[15,  16,  17].  Thus,  to  design  effective  catalytic  nano¬ 
systems  it  is  necessary  to  be  aware  of  all  the  possible 
competing  phases,  even  if  they  are  not  kinetically  ac¬ 
cessible  in  the  bulk. 

Extensive  studies  of  hafnium-based  alloys  [18],  as 
well  as  hep-based  intermetallic  alloys  [19],  demon¬ 
strated  how  high-throughput  studies  can  be  coupled 
with  other  computational  techniques,  e.g.  cluster  expan¬ 
sion,  or  with  empirical  data,  e.g.  Pettifor  maps,  to  yield  a 
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more  complete  picture  of  alloy  thermodynamics  where 
experiments  are  incomplete  or  difficult  to  carry  on. 

3.  The  Electronic  Structure  Database 

The  long  term  purpose  of  this  repository  is  to  provide 
the  community  with  electronic  structure  data  on  all  the 
known  experimental  compounds.  This  formidable  task 
is  still  in  progress.  In  the  first  stage  we  are  considering 
the  compounds  included  in  the  ICSD  database  [20] .  Of 
the  200,000+  entries  available  in  the  ICSD,  13,000  have 
already  been  calculated.  We  use  the  Aflow  materials 
framework  within  the  DFT+U  formalism,  as  explained 
in  Ref.  [8].  The  data  on  the  calculated  structures  is 
available  on  the  aflowlib.org  website  for  rapid  searching 
by  online  users.  Figure  5  shows  the  current  electronic 
structure  database  interface. 

The  user  may  specify  the  following  options  and  per¬ 
form  the  requested  search  through  the  “Submit  Query” 
button. 

•  Choose  Databases.  The  user  can  choose  the 
properties  to  be  searched  by  selecting  Structure 
[8],  Electronic  [23],  Thermoelectric  [22],  Scintil¬ 
lator  [21],  and  Magnetic  properties.  The  last  op¬ 
tion,  “Job  Status’ \  is  reserved  for  development 
purposes. 

•  Name  or  ICSD  Number.  Specify  the  composition 
of  the  compound  in  alphabetic  order,  (e.g.,  OlZnl, 
A1203)  or  the  ICSD  number  (e.g.,  82504,  9770).  If 
there  are  available  entries,  they  will  appear  in  the 
results  screen. 

•  Element(s).  Specify  a  combination  of  elements, 
using  the  keywords  &(and),  — (or),  ~(not),  '(xor), 
and  a  label  m(metals)  indicating  all  the  metals.  For 
instance,  the  user  can  specify  “Ag  &  m  &  ~Cd” 
to  obtain  all  the  multi-component  compounds  con¬ 
taining  silver,  another  metal  except  cadmium,  and 
any  other  element. 

•  Species  Number.  In  this  field  the  user  can  spec¬ 
ify  the  number  of  components  in  the  search.  For 
instance,  if  the  previous  search, 4  Ag  &  m  &  ~Cd”, 
results  in  3-,  4-,  5-  and  6-components  entries,  the 
user  can  limit  it  to  quaternary  systems  by  specify¬ 
ing  4  in  this  field. 

•  Material  Type.  In  this  field  the  user  can  choose 
to  limit  the  search  to  metals  or  semiconduc¬ 
tors/insulators. 


•  Lattice  system.  Restrict  the  results  to  one  of  the 
seven  lattice  systems:  cubic,  hexagonal,  rhombo- 
hedral,  tetragonal,  orthorhombic,  monoclinic,  tri¬ 
clinic 

•  Bravais  Lattice.  Restrict  the  results  to  one  of  the 
14  Bravais  Lattices:  cub,  fee,  bcc,  hex,  rhl,  tet,  bet, 
ore,  orcc,  oref  orci,  mcl,  mclc,  tri  [23]. 

•  Space  Group  Number.  Restrict  the  search  to 
structures  with  one  of  the  230  space  groups. 

•  Pearson  Symbol.  Specify  a  particular  Pearson 
symbol. 

•  Minimum/Maximum  band  gap.  In  these  two  fields 
the  user  can  search  for  materials  with  the  band  gap 
within  certain  limits  [21]. 

•  n-  or  p-type  Power  Factor. 

•  Minimum/Maximum  magnetic  moment.  Search 
for  materials  with  magnetic  moment/per  atom 
within  the  specified  limits. 

•  Minimum/Maximum  AS  (£/).  Search  for  materials 
with  up/down  spin  difference  at  the  Fermi  energy 
within  the  specified  limits. 

•  Aflow  version  from/to  and  Calculation  date 
from/to.  These  options  are  used  to  debug  the 
Aflow  code. 

•  Results  per  page.  The  number  of  entries  pre¬ 
sented  on  each  page  of  results. 

Typical  results  of  searches  are  shown  in  Figures  6-9. 

Figure  10  shows  an  example  of  an  entry  in  this 
database,  AliCaiOsTai  (ICSD  #99001).  The  page  in¬ 
cludes: 

•  Permanent  URL.  The  permanent  link  of  the  com¬ 
pound  on  aflowlib.org.  By  pointing  to  this  unique 
URL,  the  user  will  always  find  this  compound  in¬ 
dependently  of  the  actual  position  of  the  entry  in 
the  servers  of  the  consortium. 

•  Relaxed  Real  Space  Lattice.  Ab  initio  relaxed  ge¬ 
ometry  (< a ,  b ,  c,  a,j3,  y, ),  volume,  number  of  atoms 
per  cell  and  density  of  the  compound. 

•  Bravais  Lattice  of  the  Crystal.  Space  group,  Pear¬ 
son  Symbol,  Bravais  Lattice  of  the  primitive  cell, 
variation  of  the  Bravais  Lattice  [23],  and  Bravais 
system  for  the  crystal,  where  all  the  basis  atoms 
of  the  cell  are  considered.  This  is  the  information 
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number,  Bravais  Lattice,  Number  of  Atoms  per  Cell,  Band  Gap  (eV),  Fitted  Band  Gap  (eV),  effective  mass  of  the  electrons 
near  the  conduction  band  me,  effective  mass  of  the  electrons  at  the  minimum  of  the  conduction  band  m™in ,  effective  mass  of 
the  holes  near  the  valence  band  mh ,  effective  mass  of  the  holes  at  the  minimum  of  the  valence  band  m™in,  mismatch  mass  ratio 
mr  =  min(me/mh,mh/me ),  Valence  Band  Width  (eV),  Core  Valence  Gap  (eV),  Attenuation  Length  (cm),  Prototype  Name.  For 
definitions  see  Ref.  [21]. 


Imlrx 

tNimc 

irsn 

Bnvan 

Lattice 

Number 

of 

*Vn  >/l. 

((,W/cmK*nm) 

<P«>/L 

(jiW/cmKanm) 

<P«WL 

0iW/rmKJnm) 

<P»3>/L 

(|iW/cniK1nm) 

<Pp>/L 

(pW/cmK'nm) 

l4l 

<P|»>/L 

OiW/craK^nm) 

(41 

<Ppc>/L 

(HU/cn,K-nm) 

<Ppa>/L 

(pW/cmK'mn) 

M 

Sn 

(iiV/K) 

M 

Sp 

(uV/K) 

(4! 

Prototype  Name 

' 

F,Re(K, 

15424 

CU» 
(C-tlic  1 

5 

0.13 

0.15 

0.15 

0.15 

2.17 

2.17 

11636 

9129 

BFelKLICSD.15424 

2 

FjftiRb, 

CUB 

3 

024 

024 

0.24 

0.24 

150 

I M 

141 

141 

-91.73 

9104 

KJF€lWUUSU.49S» 

®IW 

cub 

(Qriiie.) 

3 

0JI 

031 

QJI 

Oil 

200 

100 

200 

2JOO 

■  13902 

92.92 

FelLal03.  ICSD.29il8 

Figure  8:  Results  of  a  Thermoelectric  Properties  search  within  aflowlib.org.  From  left  to  right:  index,  name,  ICSD  number,  Bravais 
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used  to  generate  the  Brillouin  Zone  [23].  Within 
Aflow  this  information  is  extracted  by  the  option 
-pgroup_xtal. 

•  Point  group  of  the  crystal.  This  information  is  im¬ 
portant  for  symmetry  characterization.  We  present 


the  crystal  family,  system  and  class,  the  point 
group  in  the  Hermann  Mauguin  and  Schonflies  no¬ 
tations,  the  orbifold,  point-group  central  symme¬ 
try,  its  order  and  its  structure.  Within  Aflow  this 
information  is  extracted  by  the  option  -edata. 
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Figure  10:  Example  of  crystallographic  and  electronic  structure  entry  in  the  aflowlib.org  database. 


•  Bravais  Lattice  of  the  Lattice.  Bravais  Lattice 
of  the  ((21,(22,(23)  lattice  vectors,  variation  of  the 
Bravais  Lattice  [23],  and  lattice  system.  This 
information  is  extracted  from  the  primitive  cells 
when  all  the  atoms  are  removed  (no-basis).  Within 
Aflow  this  information  is  extracted  by  the  option 
-pgroup. 

•  Superlattice.  If  all  the  atoms  are  considered  of  the 
same  species,  these  items  report  the  Bravais  Lat¬ 
tice,  its  variation,  the  Lattice  system  and  the  Pear¬ 
son  symbol.  This  information  is  needed  to  produce 
input  files  for  lattice  models,  such  as  the  cluster  ex¬ 


pansion. 

•  Relaxed  Reciprocal  Space  Lattice.  Geometry  of 
the  reciprocal  lattice  vectors,  volume,  Bravais  lat¬ 
tice  primitive  and  variation  [23]. 

•  Electronic  Properties.  Band  gap  (eV),  phe¬ 
nomenologically  corrected  band  gap  (eV)[21], 
band  gap  type  (direct  or  indirect),  ab  initio  total 
energy  per  atom  (eV),  magnetic  moment  per  atom 
(jUg/atom),  ab  initio  formation  enthalpy  per  atom 
(eV),  electron  mass  at  the  bottom  of  the  conduction 
band  ( me ),  hole  mass  at  the  bottom  of  the  valence 
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band  ( nth ). 

•  Band  Structure  (figure).  Band  structure  within  the 
standard  k-path  [23]. 

•  Brillouin  Zone  (figure).  Brillouin  Zone  with  the 
standard  k-path  [23]. 

•  JMOL  Structure  (figure).  A  spinning  figure  de¬ 
picting  the  atoms  in  a  unit  cell  or  supercells. 

The  information  reported  in  each  compound  entry  is 
updated  regularly  as  more  materials  properties  are  ex¬ 
tracted  from  ab  initio  calculations. 


4.  Example:  search  for  novel  nano-grained  sintered 

thermoelectric  materials 

This  database  of  electronic  properties  has  been  used 
to  search  for  improved  nano-grained  sintered  thermo¬ 
electric  materials  [22].  Thermoelectric  devices  convert 
thermal  to  electric  energy,  and  vice  versa,  with  no  mov¬ 
ing  parts.  They  are  thus  highly  reliable  and  portable 
[24],  and  are  used  in  a  variety  of  applications,  e.g.  solid- 
state  refrigerators  and  energy  harvesting  devices.  How¬ 
ever,  these  applications  are  limited  by  the  low  efficien¬ 
cies  of  the  currently  available  commercial  thermoelec¬ 
tric  compounds  [24,  25].  Amongst  the  proposed  so¬ 
lutions  to  overcome  this  limitation,  the  sintered  nano¬ 
powder  composites  are  cost-efficient  and  have  improved 
efficiency  [26].  The  grain  boundaries  in  these  materi¬ 
als  decrease  the  thermal  conductivity  k  significantly  and 
hence  enhance  the  thermoelectric  performance.  In  addi¬ 
tion  to  a  small  /c,  a  large  power  factor  P  =  cr2S ,  with  <x 
the  electric  conductivity  and  S  the  Seebeck  coefficient, 
is  required  for  high  efficiency  (characterized  by  the  di¬ 
mensionless  thermoelectric  figure  of  merit:  ZT  =  P/k , 
where  T  is  the  temperature).  The  higher  this  ZT,  the 
better  the  performance.  A  database  of  thermoelectric 
properties,  as  described  above,  was  generated  to  dis¬ 
cover  candidate  materials  a  with  large  power  factor. 

The  power  factor  was  calculated  using  first-principle 
band  structures  from  the  electronic  properties  database 
in  the  constant  mean  free  path  approximation  [22].  A 
number  of  compounds  were  found  that  may  be  promis¬ 
ing  candidates  as  nano- structured  sintered  systems.  The 
ten  n-  and  p- doped  compounds  with  highest  power  fac¬ 
tor  normalized  to  the  grain  size  in  sintered  powders  are 
shown  in  Fig.  11.  More  details  on  this  study  can  be 
found  in  Ref.  [22].  Furthermore,  a  statistical  analysis 
of  the  database  found  correlations  between  large  power 


Figure  1 1 :  The  best  ten  (a)  n-  and  (b)  p-doped  compounds  of 
highest  power  factor  normalized  by  grain  size.  The  optimal 
Seebeck  coefficients  are  shown  as  well. 


factor  and  band  gap,  effective  carrier  masses,  and  num¬ 
ber  of  atoms  in  a  unit  cell  without  any  a  priori  assump¬ 
tions.  These  correlations  can  be  used  as  guidelines  for 
finding  better  compounds  in  the  future. 

5.  Challenges  and  future  outlook.  The  Magnetic 

Database 

Currently,  the  production  of  high-performance  mag¬ 
netic  materials  depends  on  the  availability  of  few  key  el¬ 
ements.  Therefore,  the  search  for  alternative  systems  is 
one  of  the  priorities  of  the  aflowlib.org  consortium  [27]. 
We  are  developing  a  combinatorial  database  of  poten¬ 
tial  magnetic  materials.  Based  on  a  set  of  appropriate 
prototypes,  we  are  performing  an  exhaustive  scan  of  all 
multi-component  configurations.  The  results  are  then 
supplemented  with  appropriate  formation  energy  calcu¬ 
lations  to  estimate  the  feasibility  of  the  predicted  com¬ 
pounds.  Results  of  this  research  should  be  of  benefit 
to  the  magnetic  recording  industry  as  well  as  to  the  de¬ 
velopment  of  spintronics  technology.  Furthermore  our 
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Figure  12:  The  interface  of  the  magnetic  database  currently  under  development  in  aflowlib.org. 


work  may  open  new  avenues  to  the  synthesis  on  new 
high  performing  permanent  magnets  for  energy  appli¬ 
cations  [27]. 

In  Figure  12,  we  show  the  magnetic  database  search 
interface.  Similarly  to  the  electronic  structure  database 
it  offers  the  following  options: 

•  Element(s).  In  this  field  the  user  can  specify  a 
combination  of  elements.  The  syntax  is  identical 
to  the  electronic  structure  database. 

•  Prototype.  In  this  field  the  user  can  restrict  the 
search  to  particular  prototype  configurations:  all, 
AlCu2Mn,  CuL^Sn,  etc.. 

•  Band  Gap  from/to.  In  these  two  fields  the  search 
may  be  restricted  to  materials  with  band  gap  within 
the  specified  limits  [21]. 

•  Band  Gap  type.  Provides  a  choice  between  all¬ 
types,  metals,  indirect  band  gap  (I),  and  direct 
band  gap  (D). 

•  Magnetic  Moment  range  from/to.  In  these  two 
fields  the  search  can  be  restricted  to  materials  with 
a  moment/per  atom  within  the  specified  limits. 

•  Volume  range  from/to.  Restrict  the  search  to  mate¬ 
rials  with  volume/per  cell  within  the  specified  lim¬ 
its. 

•  Spin  Polarization  range  from/to.  In  these  two 
fields  the  user  can  specify  the  spin  polarization 
of  the  density  of  states  at  the  Fermi  level  (Aup  - 

N down  )  /  (^V Up  +  N^own). 


The  combination  of  database,  keywords  and  restric¬ 
tions  allow  a  flexible  search  of  potential  compounds. 
Additional  prototypes  and  options  are  planned  for  the 
near  future.  A  fully  functional,  consistent  and  robust 
magnetic  database  framework  is  expected  later  in  2012. 


6.  Summary 


We  describe  aflowlib.org,  a  repository  of  computa¬ 
tional  materials  databases  constructed  from  HT  ab  initio 
calculations  using  the  Aflow  framework. 

i.  The  metallic  binary  alloys  database  should  be  helpful 
to  the  materials  scientist  seeking  to  determine  candidate 
compounds  and  metastable  structures  in  alloy  systems 
where  such  information  is  lacking  or  controversial. 

ii.  The  electronic  structure  of  inorganic  compounds 
database  is  being  used  to  uncover  new  thermoelec¬ 
tric  and  scintillator  materials  and  should  be  helpful  in 
screening  for  additional  desirable  properties. 

iii.  The  magnetic  materials  database  would  be  used  to 
search  for  new  materials  for  magnetic,  energy  and  spin- 
tronics  applications. 

These  databases  enable  data-mining  for  trends  and 
tendencies  that  until  now  have  remained  hidden  due  to 
gaps  in  the  experimental  data  and  to  the  data  scattering 
associated  with  variability  in  the  experimental  setups 
and  methodologies.  They  should  help  overcome  some 
of  the  obstacles  associated  with  carrying  out  costly  and 
difficult  experimental  studies.  All  the  data  is  available 
the  online  repository,  aflowlib.org. 


9 


7.  Acknowledgments 

The  authors  acknowledge  Gerbrand  Ceder,  Dane 
Morgan,  Leeor  Kronik,  and  Bobby  Sumpter  for  fruitful 
discussions.  Research  supported  by  ONR  (N00014- ti¬ 
l-0136,  N00014-10-1-0436,  N00014-09-1-0921),  NSF 
(DMR-0639822,  DMR-0908753),  and  the  Department 
of  Homeland  Security  -  Domestic  Nuclear  Detection 
Office.  SC  acknowledges  support  by  the  Feinberg  fel¬ 
lowship  at  the  Weizmann  Institute  of  Science.  MBN 
wishes  to  acknowledge  partial  support  from  the  Office 
of  Basic  Energy  Sciences,  DOE  at  Oak  Ridge  National 
Laboratory  under  contract  DE-AC05-00OR22725  with 
UT-Battelle,  LLC.  SS  wishes  to  thank  Science  founda¬ 
tion  of  Ireland  (07/IN.  1/1945)  and  CRANN  for  finan¬ 
cial  support.  Computational  resources  provided  by  Ful¬ 
ton  Supercomputer  Center  at  Brigham  Young  Univer¬ 
sity,  Teragrid  (MCA-07S005),  Trinity  Center  for  High 
Performance  Computing,  Center  for  Nanophase  Mate¬ 
rials  Sciences  (CNMS20 10-206),  and  the  National  Cen¬ 
ter  for  Computational  Sciences  at  Oak  Ridge  National 
Laboratory. 

References 

[1]  P.  S.  White,  J.  R.  Rodgers,  Y.  L.  Page,  Acta  Cryst.  B  58  (2002) 
343. 

[2]  A.  Belsky,  M.  Hellenbrandt,  V.  L.  Karen,  P  Luksch,  Acta  Cryst. 
B58  (2002)  364-369. 

[3]  P.  Villars,  L.  Calvert,  Pearson’s  Handbook  of  Crystallographic 
Data  for  Intermetallic  Phases,  ASM  International,  Materials 
Park,  OH,  2nd  edition,  1991. 

[4]  P.  Villars,  K.  Cenzual,  J.  L.  C.  Daams,  F.  Hulliger,  T.  B.  Massal- 
ski,  H.  Okamoto,  K.  Osaki,  A.  Prince,  S.  Iwata,  Crystal  Impact, 
Pauling  File.  Inorganic  Materials  Database  and  Design  System , 
Binaries  Edition,  ASM  International,  Metal  Park,  OH,  2003. 

[5]  T.  B.  Massalski,  H.  Okamoto,  P.  R.  Subramanian,  L.  Kacprzak 
(Eds.),  Binary  Alloy  Phase  Diagrams,  American  Society  for 
Metals,  Materials  Park,  OH,  1990. 

[6]  P.  Villars,  M.  Berndt,  K.  Brandenburg,  K.  Cenzual,  J.  Daams, 
F.  Hulliger,  T.  Massalski,  H.  Okamoto,  K.  Osaki,  A.  Prince, 
H.  Putz,  S.  Iwata,  J.  Alloys  Compound.  367  (2004)  293-297. 

[7]  P.  Villars,  K.  Brandenburg,  M.  Berndt,  S.  LeClair,  A.  Jackson, 
Y.  H.  Pao,  B.  Igelnik,  M.  Oxley,  B.  Bakshi,  P.  Chen,  S.  Iwata,  J. 
Alloys  Compound.  317-318  (2001)  26  -  38. 

[8]  S.  Curtarolo,  W.  Setyawan,  G.  L.  W.  Hart,  M.  Jahnatek,  R.  V. 
Chepulskii,  R.  H.  Taylor,  S.  Wang,  J.  Xue,  K.  Yang,  O.  Levy, 
M.  Mehl,  H.  T.  Stokes,  D.  O.  Demchenko,  D.  Morgan,  submitted 
to  Comp.  Mat.  Sci  (2011). 

[9]  A.  Jain,  G.  Hautier,  C.  J.  Moore,  S.  P.  Ong,  C.  C.  Fischer, 
T.  Mueller,  K.  A.  Persson,  G.  Ceder,  Comp.  Mat.  Sci.  50  (2011) 
2295-2310. 

[10]  D.  de  Fontaine,  Cluster  Approach  to  Order-Disorder  Transfor¬ 
mations  in  Alloys,  volume  47  of  Solid  State  Physics,  Academic 
Press,  New  York,  1994. 

[11]  R.  H.  Taylor,  S.  Curtarolo,  G.  L.  W.  Hart,  Phys.  Rev.  B  84  (201 1) 
084101. 

[12]  O.  Levy,  R.  V.  Chepulskii,  G.  L.  W.  Hart,  S.  Curtarolo,  J.  Am. 
Chem.  Soc.  132  (2010)  833. 


[13]  O.  Levy,  M.  Jahnatek,  R.  V.  Chepulskii,  G.  L.  W.  Hart,  S.  Cur¬ 
tarolo,  J.  Am.  Chem.  Soc.  133  (2011)  158-163. 

[14]  M.  Jahnatek,  O.  Levy,  G.  L.  W.  Hart,  L.  J.  Nelson,  R.  V.  Chep¬ 
ulskii,  J.  Xue,  S.  Curtarolo,  Ordered  Structures  and  Vibrational 
Stabilization  in  Rhutenium  Alloys  from  First  Principles  Calcu¬ 
lations,  submitted  (2011). 

[15]  A.  R.  Harutyunyan,  N.  Awasthi,  A.  Jiang,  W.  Setyawan, 
E.  Mora,  T.  Tokune,  K.  Bolton,  S.  Curtarolo,  Phys.  Rev.  Lett. 
100  (2008)  195502. 

[16]  S.  Curtarolo,  N.  Awasthi,  W.  Setyawan,  A.  Jiang,  K.  Bolton, 
T.  Tokune,  A.  R.  Harutyunyan,  Phys.  Rev.  B  78  (2008)  054105. 

[17]  R.  V.  Chepulskii,  S.  Curtarolo,  ACS  Nano  5  (2011)  247-254. 

[18]  O.  Levy,  G.  L.  W.  Hart,  S.  Curtarolo,  Acta  Mater.  58  (2010) 
2887-2897. 

[19]  O.  Levy,  G.  L.  W.  Hart,  S.  Curtarolo,  Phys.  Rev.  B  81  (2010) 
174106. 

[20]  A.  D.  Mighell,  V.  L.  Karen,  Acta  Cryst.  A49  (1993)  c409. 

[21]  W.  Setyawan,  R.  M.  Gaume,  S.  Lam,  R.  S.  Feigelson,  S.  Cur¬ 
tarolo,  ACS  Comb.  Sci.  13  (201 1)  382-390. 

[22]  S.  Wang,  Z.  Wang,  W.  Setyawan,  N.  Mingo,  S.  Curtarolo,  Phys. 
Rev.  X  (2011). 

[23]  W.  Setyawan,  S.  Curtarolo,  Comp.  Mat.  Sci.  49  (2010)  299. 

[24]  G.  S.  Nolas,  J.  Sharp,  H.  J.  Goldsmid,  Thermoelectrics:  Basic 
Principles  and  New  Materials  Developments,  Springer  Verlag, 
2001. 

[25]  C.  B.  Vining,  Nat.  Mater.  8  (2009)  83-85. 

[26]  G.  Joshi,  H.  Lee,  Y.  Lan,  X.  Wang,  G.  Zhu,  D.  Wang,  R.  W. 
Gould,  D.  C.  Cuff,  M.  Y.  Tang,  M.  S.  Dresselhaus,  G.  Chen, 
Z.  Ren,  Nano  Lett.  8  (2008)  4670-4674. 

[27]  N.  Jones,  Nature  472  (2011)  22-23. 


10 


